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Alternative Scenarios of Relativistic Heavy-Ion Collisions: 
II. Particle Production 

Yu.B. Ivanoyi-El 
' Kurchatov Institute, Moscow RU-123182, Russia 

Particle production in relativistic collisions of heavy nuclei is analyzed in a wide range of incident 
energies 2.7 GeV < y^sjviv < 62.4 GeV. The analysis is performed within the three-fluid model 
employing three different equations of state (EoS): a purely hadronic EoS, an EoS with the first- 
order phase transition and that with a smooth crossover transition. It is found that the hadronic 
scenario fails to reproduce experimental yields of antibaryons (strange and nonstrange), starting 
already from lower SPS energies, i.e. ^/Jnn > 5 GeV. Moreover, at energies above the top SPS one, 
i.e. ^/snn > VTA GeV, the mid-rapidity densities predicted by the hadronic scenario considerably 
exceed the available RHIC data on all species. At the same time the deconfinement-transition 
scenarios reasonably agree (to a various extent) with all the data. The present analysis demonstrates 
certain advantage of the deconflnement-transition EoS's. However, all scenarios fail to reproduce 
the strangeness enhancement in the incident energy range near 30^4 GeV (i.e. a horn anomaly in 
the _ft'+/7r+ ratio) and yields of (^-mesons at 20A-40^ GeV. 

PACS numbers: 25.75.-q, 25.75.Nq, 24.10.Nz 
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I. INTRODUCTION 

This paper continues a series of reports on simulations 
of relativistic heavy-ion collisions within different scenar- 
ios [l|, |2]- These simulations were performed within a 
model of the three-fluid dynamics (3FD) Q employing 
three different equations of state (EoS): a purely hadronic 
EoS 3 (hadr. EoS), which was used in the major part 
of the 3FD simulations so far [3, M^ , and two versions 
of EoS involving the deconfinement transition Q . These 
two versions are an EoS with the first-order phase transi- 
tion and that with a smooth crossover transition. Details 
of these calculations are described in the first paper of 
this series Q dedicated to analysis of the baryon stop- 
ping. 

In this paper I report results on particle production 
and rapidity distributions of these particles in relativis- 
tic heavy-ion collisions in the energy range from 2.7 
GeV to 62.4 GeV^ in terms of the center-of-mass in- 
cident energy {^/s^n)- This domain covers the en- 
ergy range of the RHIC (Relativistic Heavy- Ion Collider) 
beam-energy scan and SPS (Super Proton Synchrotron) 
low-energy-scan programs, as well as energies of the fu- 
ture FAIR (Facility for Antiproton and Ion Research) 
and NICA (Nuclotron-based Ion Collider Facility) facil- 
ities and the AGS (Alternating Gradient Synchrotron) 
at BNL (Brookhaven National Laboratory). As demon- 
strated in the first papers of this series [l|, y] , within the 
considered here first-order-transition and crossover sce- 
narios the deconfinement transition takes place in the 
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^ Results for the top calculated energy of 62.4 GeV should be taken 
with care, because they are not quite accurate, as an accurate 
computation requires unreasonably high memory and CPU time. 



region of top-AGS-low-SPS incident energies. The ex- 
perimental baryon stopping also indicates certain signs 
of a deconfinement transition. Therefore, in this paper 
the attention is primaraly focused on this incident en- 
ergy range. It should be mentioned that available data 
on particle production in the AGS-SPS energy range have 
already been analyzed within various models |3|, |5|, llOl - 
llTJ I . The aim of this paper is a comparative analysis of 
these data within different assumptions on the EoS. 

Closely related to the particle production, the hadron 
yield ratios are also discussed. The most intriguing issue 
among them is a horn anomaly in the K~^ /tt^ ratio |37l . 
ISSf that still has no explanation within any dynamical 
model. Unfortunately, the 3FD model is not an exception 
from this list. 



II. 3FD MODEL 

The 3FD model @ is a straightforward extension of 
the 2-fluid model with radiation of direct pions J18l - [20l | 
and (2-|-l)-fluid model [2l|, [22l. The above models were 
extend in such a way that the created baryon-free fluid 
(which is called a "fireball" fluid, following the Frankfurt 
group) is treated on equal footing with the baryon-rich 
ones. A certain formation time r is allowed for the fire- 
ball fluid, during which the matter of the fluid propagates 
without interactions. The formation time is associated 
with a flnite time of string formation. It is similarly in- 
corporated in kinetic transport models such as UrQMD 
m andHSD JH. 

Unlike the conventional hydrodynamics, where local 
instantaneous stopping of projectile and target matter is 
assumed, a specific feature of the 3FD is a finite stopping 
power resulting in a counter-streaming regime of leading 
baryon-rich matter. The basic idea of a 3-fluid approxi- 
mation [1^, [25| is that at each space-time point a gener- 



ally noncquilibrium distribution function of baryon-rich 
matter can be represented as a sum of two distinct con- 
tributions initially associated with constituent nucleons 
of the projectile (p) and target (t) nuclei. In addition, 
newly produced particles, populating the mid-rapidity re- 
gion, are associated with a fireball (f) fiuid. Therefore, 
the 3-fluid approximation is a minimal way to simulate 
the finite stopping power at high incident energies. Each 
of these fluids is governed by conventional hydrodynamic 
equations which contain interaction terms in their right- 
hand sides. These interaction terms describe mutual fric- 
tion of the fluids and production of the flreball fluid. The 
friction between fluids was fitted to reproduce the stop- 
ping power observed in proton rapidity distributions for 
each EoS, as it is described in Ref. [21 in detail. 

A conventional way of applying the fluid dynamics to 
heavy-ion collisions at RHIC and LHC energies is to pre- 
pare the initial state for the hydrodynamics by means of 
various kinetic codes [111, [l3, ll3, 123 • Contrary to these 
approaches, the 3FD model treats the collision process 
from the very beginning, i.e. the stage of cold nuclei, 
up to freeze-out within the fluid dynamics. The con- 
ventional hydrodynamical model of Refs. [23, UM does 
such simulations in a very similar way but without tak- 
ing into account incomplete stopping of colliding nuclei 
at the initial stage of the reaction. Therefore, such kind 
of simulations are justified only at moderately high ener- 
gies. 

Freeze-out is performed accordingly to the procedure 
described in Ref. 3 and in more detail in Refs. [29l[30|. 
This is a modified Milckhin version of the freeze-out 
that possesses exact conservation of the energy, mo- 
mentum and baryon number. Contrary to the conven- 
tional Cooper-Frye approach [3l| , the modified Milekhin 
method has no problem associated with negative con- 
tributions to particle spectra. The freeze-out criterion 
is based on a phenomenologically determined freeze-out 
energy density etrz- This method of freeze-out can be 
called dynamical, since the freeze-out process here is inte- 
grated into fiuid dynamics through hydrodynamic equa- 
tions. This kind of freeze-out is similar to the model of 
"continuous emission" proposed in Ref. JS^]. There the 
particle emission occurs from a surface layer of the mean- 
free-path width. In the 3FD case the physical pattern is 
similar, only the mean free path is shrunk to zero. To 
the moment of the freeze-out the matter is already in the 
hadronic phase in the case of the 2-phase EoS, while for 
the crossover EoS this is not so. However, this is not 
a problem because, any case, the thermodynamic quan- 
tities of the frozen-out matter are recalculated from the 
in-matter EoS, with which the hydrodynamic calculation 
runs, to the hadronic gas EoS. This is done because a part 
of the energy is still accumulated in collective mean fields 
at the freeze-out instant. This mean-field energy should 
be released before calculating observables. Otherwise, 
the energy conservation would be violated. 

The baryon stopping turns out to be only moderately 
sensitive to the freeze-out energy density. The freeze- 



out energy density eirz = 0.4 fm/fm^ was chosen mostly 
on the condition of the best reproduction of secondary 
particles yields. 

In the 3FD model ^ , particles are not isotopically dis- 
tinguished, i.e. the model deals with nucleons, pious, etc. 
rather than with protons, neutrons, tt^, tt^, tt^, etc. In 
fact, it is not a problem to formulate the model for iso- 
topically distinguished species. For that it is necessary 
to introduce into the EoS an additional isotopic chemi- 
cal potential associated with the electric-charge conser- 
vation, and an additional continuity equation for each 
fiuid that controls this electric-charge conservation. The 
problem comes out at the stage of numerical solution of 
this set of equations. This is the curse of dimensionality. 
Treatment of a 3-dimensional EoS (temperature and two 
chemical potentials) together with the (3-f l)-space-time 
evolution is too complicated for computers available at 
present. Thus, the disregard of the isospin is a forced 
approximation. 

Presently it is unconventional to apply hydrodynamics 
to such low incident energies as 2A GeV. However, 20-30 
years ago the hydrodynamics was a conventional tool for 
analysis of the heavy-ion collisions at incident energies 
around lA GeV, see e.g. Refs. [H, [H, [H (H . Ah the 
arguments put forward in favor of the hydrodynamic 
approach in those references still hold true. In particular, 
the popular now blast- wave model [35l . |36| was originally 
developed for these low energies. 

III. RAPIDITY DISTRIBUTIONS OF 
PRODUCED PARTICLES 



The rapidity distributions of produced particles were 
calculated at fixed impact parameters (6). Correspon- 
dence between experimental centrality, i.e the fraction of 
the total reaction cross section related to a data set, and 
a mean value of the impact parameter was taken from 
the paper [SOj in case of NA49 data. For Au+Au col- 
lisions it was approximately estimated proceeding from 
geometrical considerations. Feed-down from weak decays 
of strange particles into non-strange hadron yields was 
disregarded in accordance with measurement procedure 
of the NA49 collaboration. At the AGS energies this con- 
tribution of weak decays is just negligible. At RHIC en- 
ergies (Au-|-Au collisions at beam-energy-scan energies), 
the contribution of weak decays into non-strange hadron 
yields was taken into account similarly to that done in 
data of the STAR and PHENIX coUaborations. 

In Fig. [U comparison with available experimental data 
on rapidity distributions of pious from central collisions is 
presented. As particles are not isotopically distinguished 
in the 3FD model [3| , only the total number of pious Nt^ 
is calculated in the model. In Fig. [1] rapidity distribu- 
tion of Nt^/S pious is displayed which, strictly speaking, 
should be compared with the sum of 7r+ , 7r~ and 7r° dis- 
tributions divided by 3. In view of the fact that the 
experimental tt~ distributions are always situated above 
the TT"*" ones, a "good agreement" with the model means 
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FIG. 1: (Color online) Rapidity spectra of pions from central collisions of Au+Au {b — 2 fm) and Pb+Pb (6 = 2.4 fm). 
Experimental data are taken for central Au+Au collisions at AGS energies [401 ] (at 2A-8A GeV, centrality 5%) and [4l[| (at 
10.7 A GeV, centrality 3%), and Pb+Pb collisions at SPS energies [H.liil with centrality 7% at 20^-80^ GeV and centrality 
5% at 15SA GeV. 



that calculated spectrum of 1/3 of all pions turns out in 
between experimental n~^ and tt~ distributions. As seen 
from Fig. [1] this is always the case at AGS energies for 
all three EoS's. At the SPS energies differences between 
predictions of different scenarios and experimental data 
do not exceed 10%, which is quite acceptable in view of 
uncertainties of the model (such as a choice of the impact 
parameter, accuracy of the multi-fluid approximation, 
phenomenological freeze-out procedure, etc.). Therefore, 
we can approximately conclude that the pion rapidity 
distributions are reasonably well reproduced by all of the 
considered EoS's in the region of AGS-SPS incident en- 
ergies. 

A comment concerning the hadronic EoS is in order 
here. As it was mentioned in Refs. [I, Q, we pro- 
ceeded from a principle of fair treatment of any EoS. 
It means that any possible uncertainties in the parame- 
ters are treated in favor of the EOS. Therefore, for the 
hadronic EoS the parameters of the friction enhancement 
and the formation time of the fireball fluid l2,y| were cho- 
sen such that the baryon stopping and the mid-rapidity 
pion density were reproduced at the top SPS energy of 
158^4 GeV. Thus, it is not surprising that the hadronic- 
EoS result is the best at 158A GeV, see Fig. [TJ However, 
even with this fit the pion distributions at midrapidity, 
i.e. (7r+ + 7r^)/2; turn out to be somewhat underesti- 
mated at lower SPS energies, as seen in Fig. [TJ and 
strongly overestimated at higher (RHIC) energies, as it 



will be seen below. 

In Fig. [51 rapidity distributions of positive and nega- 
tive kaons from central collisions are presented. In fact, 
the model provides distributions of kaons, K, and an- 
tikaons, K, which are sums of yields of K'^ + K'^ and 
K^ + K'^, respectively. The distributions of K'^ and K^ 
are then approximately associated with distributions of 
K and K, respectively, divided by two. This approxi- 
mation induces a certain error especially at low incident 
energies, where yields of K'^ and K'^ (and, correspond- 
ingly, K~ and K'^) noticeably differ. 

As seen from Fig. [21 the K^ production is distinctly 
overestimated in all of the considered scenarios. How- 
ever, it is early to conclude that all the scenarios fail. 
This overestimation can be a consequence of the fact 
that the extrapolation of experimental data (the expo- 
nential one) beyond the measured points (thin solid lines 
in Fig. [3]) essentially differs from results of calculations 
at high transverse masses (jut), while in the experimen- 
tally measured Iow-to^ region these are quite similar. In 
fact, determination of the rapidity density requires inte- 
gration over TTiT up to infinity. Therefore, the experimen- 
tal result can be different from the model ones because 
of their difference in the mx region, where experimental 
data are absent. Similar situation takes place for the if "*" 
yields. At lower SPS energies the experimental extrap- 
olation to high TTiT reasonably agree with model results 
and hence the same reasonable agreement is observed in 




FIG. 2: (Color online) Rapidity spectra of positive [(a) block of panels] and negative [(b) block panels] kaons from central 
collisions of Au+Au (6=2 fm) and Pb+Pb (6 = 2.4 fm). Exp erimental data are for central Au+Au collisions at AGS energy 
of 10.74 GeV (centrality 3%) ^ and Pb+Pb collisions [sl El] at SPS energies with centrality 7% at 20A-80A GeV and 
centrality 5% at 1584 GeV. 
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FIG. 3: (Color online) Transverse mass spectra at mid-rapidity of positive (a) and negative (b) kaons from central collisions 
Pb+Pb (6 = 2.4 fm) at various incident energies and those in the extended transverse-mass range at 158^4 GeV (c). Experimental 
data are taken from Refs. [38 . |4^ | and for high pr momenta, from Ref. [43] . 
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FIG. 4: (Color online) Rapidity spectra of (j> mesons from central collisions of Au-|-Au (6 = 2 fm) and Pb-|-Pb {b = 2.4 fm). 
Experimental data for central Pb-|-Pb collisions at SPS energies with centrality 7% at 20A-80A GeV and centrality 5% at 158A 
GeV are taken Ref. [451 ] and for central Au+Au collisions at 11.7 A GeV/c and centrality 5%, from Ref. [4a]. 
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FIG. 5: (Color online) Rapidity spectra of A [(a) block of panels] and A [(b) block of panels] hyperons from central collisions of 
Au+Au (6 = 2 fm) and Pb+Pb (b = 2.4 fm). Results for Pb+Pb collisions with b = 4.6 fm are also displayed for A's at 158 A 
GeV. Experimental data for central Pb+Pb collisions at SPS energies with centrality 7% at 20^4-80^ GeV and centrality 10% 
at 158A GeV are taken from Ref. [4J], and for central Au+Au collisions at 11.5^ GeV/c and centrality 4%, from Ref. [471 ]. 
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FIG. 6: (Color online) Rapidity spectra of S [(a) block of panels] and H+ [(b) block of panels] hyperons from central collisions 
of Au+Au (& = 2 fm) and Pb+Pb {b = 2.4 fm). Experimental data are taken for central Pb+Pb collisions [4J] at SPS energies 
with centrality 7% at 20A-80A GeV and centrality 10% at 158A GeV. Results for Pb+Pb collisions with b = 4.6 fm are also 
displayed for H^'s at 158A GeV. 




FIG. 7: (Color online) Rapidity spectra of anti-protons from central collisions of Au+Au (6=2 fm) and Pb+Pb (6 = 2.4 fm). 
Experimental data for central Pb+Pb collisions at 158 A GeV with centrality 5% are taken from Ref. pa ]. 
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FIG. 8: (Color online) Rapidity spectra of anti-protons at 
various centralities from collisions Pb-fPb at _B;af, = 158^1 
GeV. Experimental data are taken from Ref. [4g|. Curves 
present calculations with the crossover EoS at different impact 
parameters (6). 

rapidity distributions. At higher SPS energies (80A and 
158A GeV) the experimental (exponential) extrapola- 
tion is again distinctly lower than the model predictions. 
Therefore, the model rapidity distributions overestimate 



the experimental ones. It is worthwhile to notice that the 
calculated shapes of the ttit spectra of kaons are some- 
what similar to that of pion spectra, i.e. they are slightly 
concave in the logarithmic scale rather than linear, as 
it was assumed in the experimental extrapolation. Mea- 
surements of high transverse momenta later performed 
by NA49 collaboration at 158 A GeV ^, showed that 
the exponential extrapolation of the low-mr data indeed 
underestimates the high-TTiT data, see (c) panel of Fig. 
[31 At the same time the 3FD predictions overestimate 
these high-TOT data. The latter fact is expected, because 
even abundant hadronic species become rare species at 
high momenta. Therefore, their treatment on the basis of 
grand canonical ensemble (which is the case in the 3FD 
model) results in overestimation of their yield. 

In view of the above possible explanation of the differ- 
ence between calculated and experimental spectra it is 
early to conclude that the 3FD model fails to reproduce 
the K'^ and K~ rapidity distributions. Again different 
EoS-scenarios agree (or disagree) with the available data 
approximately to the same extent. An important feature 
of the K^ distributions is that all scenarios underesti- 
mate mid-rapidity values of the experimental data. This 
results in a failure to reproduce the experimental "horn" 
in the if+/7r+ ratio [33l, see sect. IIVI 

In Fig. [5] rapidity distributions of mesons from cen- 
tral collisions are presented. Here the situation is contro- 
vcrsal. At high SPS energies, predictions of the hadronic 
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FIG. 9: (Color online) Rapidity densities dN/dy at mid-rapidity of various produced particles as functions of the center-of-mass 
energy of colliding nuclei predicted by 3FD calculations with three considered EoS's. Experimental data are from compilation 
of Ref. [49II complemented by recent data from STAR collaboration [50] and latest update of the compilation of NA49 numerical 
results [5ll.[52l. 



EoS look certainly preferable. At lower SPS energies and 
the top AGS energy, predictions of different scenarios 
are rather similar and all of them considerably underes- 
timate the experimental data at the mid-rapidity region. 
A reason of this is still unclear. 

Let us turn to baryonic distributions. Since proton dis- 
tributions have been already analyzed in detail in Refs. 
[lljQl, let us directly proceed to strange baryons. In Figs. 
[5] and [S] rapidity distributions of A, A, S~ and S"*" hyper- 
ons from central collisions are presented. Because of the 
low centrahty at 158A GeV (10%), results for Pb+Pb 
collisions with b = 4.6 fm are also displayed for A and 
S^. For corresponding antiparticles A and S+ the re- 
sults with h = 2.4 fm and b = 4.6 fm are fairly close to 
each other. Therefore, calculations with b — 4.6 fm are 
not presented. As seen, data on A and S^ hyperons are 
reasonably reproduced by all of considered EoS's. Ap- 
parent preference of the hadronic scenario in reproduc- 
ing the S~ data cannot be considered as an argument in 
favor of this scenario because the extent of the overesti- 
mation by the deconfinement-transition scenarios is just 
within the accuracy of the grand-canonical approach. At 
the same time, for anti-hyperons A and S+ the hadronic 



scenario evidently fails: it considerably overestimates the 
data, especially at the top SPS energy. The case of S+ 
at 30^ GeV is not spectacular. There the data have too 
large error bars. For A and A the crossover EoS looks 
certainly preferable. At the same time it is hardly pos- 
sible to expect from the 3FD model with 2-phase and 
crossover EoS's a better agreement with data on such 
rare probes as S~ and S"*" because the model is based on 
grand-canonical statistics and therefore requires "large" 
multiplicities to be valid. 

Figure [7] confirms that the hadronic scenario strongly 
overestimates experimental data on antibaryons, in the 
present case, on antiprotons. At the same time, the 
crossover scenario almost perfectly reproduces the an- 
tiproton rapidity distribution even at various centrali- 
ties, see Fig. [51 Correspondence between the centrahty 
related to a data set and a mean value of the impact pa- 
rameter used in the calculation is taken from Ref. [30] . 

In all the blocks of figures I kept the results for Au-f Au 
collisions at 10 A GeV, even if there are no data for this 
reaction. This done because the deconfinement transi- 
tion in the presently studied scenarios of nuclear colli- 
sions starts around this incident energy, see Refs. [l|, [2|. 
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Therefore, it is important to observe evolution of particle 
distributions beginning from this energy. 
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FIG. 10: (Color online) Strangeness suppression factor as a 
function of the center-of-mass energy of colliding nuclei. 



IV. MID-RAPIDITY DENSITIES AND THEIR 
RATIOS 

Fig. [9] summarizes and extends the discussion of the 
previous section. It displays excitation functions of the 
mid-rapidity values of the rapidity spectrum of various 
produced particles but in a wider range of incident en- 
ergies Y^J/vjv = 2.7-62.4 GeV. Results for the top cal- 
culated energy of ^/snn = 62.4 GeV should be taken 
with care, because they are not quite accurate, as an ac- 
curate computation requires unreasonably high memory 
and CPU time. 

The strangeness production at low incident energies is 
overestimated within the 3FD model. This is the con- 
sequence of the fact that the 3FD model is based on 
the grand canonical ensemble. This shortcoming can 
be easily cured by introduction of a phcnomenological 
factor 75 5J] which takes into account an additional 
strangeness suppression due to constraints of canonical 
ensemble. The mid-rapidity densities of single-strange 
particles {K^, A and A), displayed in Fig. [9l are mul- 
tiplied by 7s factor, and multi-strange particles {(j), S", 
S"*" and fi), by (75)"" factor, where Ug is the number of 
valence strange quarks contained in this particle. 

The excitation function of the 75 factor is presented in 
Fig. [TUl which is of course applicable only to central col- 
lisions of considered nuclei. On average, there is no need 
for additional strangeness suppression at y/sNN > 5 GeV, 
though the reproduction of data on some strange species 
is still far from being perfect. The hadronic scenario looks 
worst in this respect. At 5 GeV < y^SNN < 17 GeV, i.e. 



in the range of the SPS energy, it is impossible to find a 
unique value of 75 for this scenario to simultaneously im- 
prove reproduction of data on all of the strange hadrons. 
At ^ysNN > 17 GeV the hadronic scenario again requires 
an additional strangeness suppression. However, intro- 
duction of such a suppression at high incident energies 
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FIG. 11: (Color online) The energy dependence of ratios of 
the rapidity densities dN/dy at mid-rapidity: (a) K'^ to tt^, 
(b) K^ to n~ , and (c) A to 7r~. Compilation of experimental 
data is from Ref. [52| complemented by recent STAR data 
[sol . [53 . Calculated ratios with pions including contribution 
from weak decays are also displayed by thin lines. The latter 
ratios are relevant to STAR and PHENIX data. 



is physically unjustified, because multiplicities of strange 
hadrons are already quite high for applicability of the 
grand canonical ensemble. 

Fig. [9] confirms that the hadronic scenario consider- 
ably overestimates mid-rapidity densities of antibaryons 
(strange and nonstrange) predicted by EoS's with a dc- 
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confinement transition and (as a rule) those deduced form 
experiment, starting already from lower SPS energies, 
i.e. x20A GeV. At energies above the top SPS one, the 
mid-rapidity densities predicted by the hadronic scenario 
considerably exceed the data on all species. Notice that 
agreement of hadronic-scenario predictions with data on 
the major part of species in the SPS energy range was 
achieved at the expense of considerable enhancement the 
inter-fluid friction in the hadronic phase 0, 3 a s com- 
pared with its microscopic estimate of Ref. [55[. Pre- 
cisely this enhancement makes bad job at energies above 
the SPS. 

The advantage of deconfinement-transition scenarios is 
that they do not require any modification of the micro- 
scopic friction in the hadronic phase. They reasonably 
agree with all the data with the exception of negative 
kaons and 0-mesons in the SPS energy range. A possible 
reason of disagreement with negative-kaon data is inaccu- 
racy of the experimental extrapolation of the transverse- 
mass spectra, discussed in the previous section. As for 
0-mesons, the reason of this is still unclear. Rare probes, 
such as S"'" and f2~, are also poorly reproduced, that is 
a consequence of inapplicability of the grand canonical 
ensemble for the treatment of such rare species. 

In Fig. Ill[ ratios of various hadronic mid-rapidity den- 
sities are presented as functions of the incident energy. 
Calculated ratios with pions including contribution from 
weak decays of strange hadrons are also displayed by the 
corresponding thin lines. These ratios are relevant to 
STAR and PHENIX data [13, IH. All of the consid- 
ered scenarios do not reproduce the horn anomaly in the 
K+/TT+ ratio ^ around Eiab = 30 A GeV. As it was 
demonstrated in Fig. [2l the reason is that all scenarios 
underestimate mid-rapidity values of K'^ as compared 
to the experimental data, while its rapidity distributions 
at peripheral rapidities look almost perfect. As it was 
found, the value of this ratio is quite insensitive to varia- 
tions of the freeze-out parameter Sfrz 0, 9 , because the 
yields of K~^ and tt"*" increase or decrease proportionally 
with the change of £{rz ■ 

Authors of Ref. [53 have noticed that within the sta- 
tistical model the peak the K~^ /tt'^ ratio can be better 
reproduced if the contribution of the /o(600) resonance, 
i.e. the a meson, is taken into account. However, au- 
thors of Ref. [531 concluded that the a meson does not 
solve the problem. In the present 3FD calculations the 
a meson is included. Moreover, the a spectral function 
is taken into account accordingly to parametrization of 
Ref. m. 

A possible way to improve the reproduction of the horn 
anomaly is to accept different freeze-out conditions, i.e. 
energy densities Sfrz, for K~^ and tt"*" in this narrow^ en- 
ergy region near Eiab = 30A GeV. However, physical 



^ In a wide energy range, the universal freeze-out for K^ and tt^ 
works faily well and there are no reasons to change it. 



reasons for such a radical solution of the horn-anomaly 
problem are absent. Another possibility consists in no- 
ticeable difference in K~^ and K^ densities in this energy 
range. I would like to remind that the 3FD model does 
not distinguish isotopic states of hadrons and hence the 
yield of K^ is determined as half of the total K yield 
consisting of the sum of K^ and K'^ contributions. Simi- 
larly, the 7r+ yield is determined as one third of the total 
pion one. However, a physical reason for such an iso- 
topic asymmetry in this narrow energy range is not seen. 
Moreover, such an asymmetry would be opposite to a 
natural isotopic asymmetry at which the K'^ yield ex- 
ceeds the K^ one. This natural asymmetry results from 
initial isotopic asymmetry of colliding nuclei. 

A possible reason for overcstimation of K~ /■k~ ratio, 
see (b) panel of Fig. [TTl as compared with the data has 
been already mentioned in the previous section. It can 
result from difference of the extrapolation of transverse- 
mass experimental data (exponential extrapolation) be- 
yond the measured points (thin solid lines in Fig. ^ from 
results of calculations at high transverse masses (mx). 
Therefore, it would be desirable to directly compare the 
transverse-momentum spectra predicted by the model 
with experimental ones at high transverse momenta. The 
A/tt^ ratio, see (c) panel of Fig. [TTl is quite satisfactorily 
reproduced by all of the scenarios. 



V. SUMMARY 

Results on rapidity distributions of particles produced 
in relativistic heavy-ion collisions in the energy range 2.7 
GeV < y/SNN < 62.4 GeV are presented. These simu- 
lations were performed within the three-fluid model [3| 
employing three different equations of state: a purely 
hadronic EoS [Jl and two versions of the EoS involving 
the deconfinement transition [9[. These two versions are 
an EoS with the first-order phase transition and that with 
a smooth crossover transition. Details of these calcula- 
tions are described in the first paper of this series [2| 
dedicated to analysis of the baryon stopping. 

If was found that the model within all of the consid- 
ered scenarios reproduces approximately to the same ex- 
tent the available data on almost all (with some excep- 
tions) hadronic species in the AGS-SPS energy range, 
i.e. at 2.7 GeV < y/sWw < 17.4 GeV. In the case of 
the hadronic EoS this agreement is achieved at the ex- 
pense of noticeable enhancement of the inter-fluid friction 
in the hadronic phase [2|, [3i as compared with its micro- 
scopic estimate of Ref. [53 ■ Reproduction of the data is 
better for abundant species, as it should be within grand- 
canonical statistics. The above mentioned exceptions are 

• At lower SPS energies and the top AGS energy, pre- 
dictions of all scenarios are rather similar and all of 
them considerably underestimate the experimental 
0-meson data. Moreover, at the top SPS energy, 
the hadronic EoS unexpectedly gives the best de- 
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scription of the 
is unclear. 



-meson data. The reason of this 



• The strangeness enhancement in the incident en- 
ergy range near 30A GeV, i.e. the "horn" anomaly 
in the K~^/t:^ ratio |38l |. is not reproduced in any 
scenario. 

• The K~ production is distinctly overestimated in 
all of the considered scenarios. This overestimation 
can be a consequence of the extrapolation of ex- 
perimental data (the exponential one) beyond the 
measured points that essentially differs from results 
of calculations at high transverse masses, while in 
the experimentally measured Iow-ttit^ region these 
are quite similar. 

• Data of rare probes (like S^ and S+ hyperons) are 
poorly reproduced apparently because of the grand- 
canonical statistics used in the 3FD model. 

• The last exception concerns only the hadronic 
EoS. The hadronic scenario considerably over- 
estimates experimental rapidity densities of an- 
tibaryons (both strange and nonstrange), starting 
already from lower SPS energies, i.e. « 20A GeV. 
At the same time, the deconfinement-transition sce- 
narios successfully (to a various extent) reproduce 
these data. 

The failure with the t/f-meson data and the "horn" 
anomaly requires further investigation. The proposed 
reason of the overestimation of the K^ production also 
needs to be checked. 

At energies above the top SPS one, i.e. at ^/snn > 
17.4 GeV, the mid-rapidity densities predicted by the 
hadronic scenario considerably exceed the available 
RHIC data on all the species. The deconfinement- 
transition scenarios reasonably agree with all the data. 

The 3FD model deals with three different fluids [l,!]. 
These are two baryon-rich fluids initially associated with 
constituent nucleons of the projectile and target nuclei. 
These fluids are either spatially separated or unified at 
the freeze-out stage. In addition, newly produced parti- 
cles, populating the mid-rapidity region, are associated 
with a baryon-free ( "fireball" ) fluid which remains undis- 
solved in baryonic fluids till the freeze-out. A certain 
formation time r is allowed for the fireball fluid, during 
which the matter of the fluid propagates without inter- 
actions. The formation time is associated with a finite 
time of string formation. The main difference concerning 
this baryon-free fiuid in considered alternative scenarios 
consists in different formation times: r = 2 fm/c for the 
hadronic scenario and t = 0.17 fm/c for scenarios involv- 
ing the deconfinement transition |2|. 

As it is seen from simulations, the main contribution 
to baryon and meson yields comes from baryon-rich fiu- 
ids within the deconfinement-transition scenarios. This 



contribution slowly decreases with the incident energy 
rise. For instance, at y^SNN = 39 GeV approximately 
half of pious are produced from the baryon-free fluid 
within the deconfinement-transition scenarios. For all 
other particles (i.e. baryons and mesons) and lower en- 
ergies this fraction is essentially lower. At the same time 
the fraction of half for pious from the baryon-free fluid 
is achieved already at -y/s]v7V — 9 GeV (i.e. Eiab — 40A 
GeV) within the hadronic scenario. This is one of the 
reasons why r was chosen so large in the hadronic sce- 
nario. Large formation time prevents absorption of the 
baryon-free matter by the baryonic fluids. Without this 
large contribution of the baryon-free fluid it is impossi- 
ble to reproduce mesonic yields at SPS energies within 
the hadronic scenario. However, this strongly developed 
baryon-free fluid makes bad job for antibaryons in the 
case of the hadronic EoS. The reason is that antibaryons 
are dominantly produced from the baryon-free fluid even 
at lower incident energies within all of the considered 
scenarios. Their yields in the hadronic scenario strongly 
overestimate experimental data. At y^iivlv > 17.4 GeV 
the overdeveloped baryon-free fluid makes already bad 
job for for all produced particles, resulting in consider- 
able overestimation of the available RHIC data on all the 
species within the hadronic scenario. 

All the above speculations are formulated in terms of 
the 3FD model. However, they hold true in general phys- 
ical terms, if baryon-rich fluids are understood as a mat- 
ter formed of leading particles, predominantly (but not 
solely) occupying peripheral rapidities, and the baryon- 
free fluid is associated with newly produced particles, 
predominantly (but not solely) occupying mid-rapidity 
region. At large baryon stopping, the difference between 
these groups of particles become less pronounced. 

All the above results certainly indicate in favor of the 
deconfinement-transition scenarios, though they also are 
not perfect in reproduction of the data ((/)-meson, the 
"horn" anomaly). It is difficult to judge on a prefer- 
ence of the 2-phase or crossover EoS's. The first-order 
transition in the 2-phase EoS is very abrupt. On the 
contrary, the crossover transition constructed in Ref. [9| 
is very smooth [l|, Q. In this respect, this version of the 
crossover EoS certainly contradicts results of the lattice 
QCD calculations, where a fast crossover, at least at zero 
chemical potential, was found [53 . Therefore, a true EoS 
is somewhere in between the crossover and 2-phase EoS's 
of Ref. 0. 
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